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Abstract — We investigate the secure connectivity of wire- 
less sensor networks under the pairwise key distribution 
scheme of Chan et al.. UnUke recent work which was 
carried out under the assumption of full visibility, here we 
assume a (simplified) communication model where unreli- 
able wireless hnks are represented as on/off channels. We 
present conditions on how to scale the model parameters so 
that the network i) has no secure node which is isolated and 
ii) is securely connected, both with high probability when 
the number of sensor nodes becomes large. The results are 
given in the form of zero-one laws, and exhibit significant 
differences with corresponding results in the full visibiUty 
case. Through simulations these zero-one laws are shown to 
be valid also under a more reaUstic communication model, 
i.e., the disk model. 

Keywords: Wireless sensor networks. Security, Key predis- 
tribution. Random graphs. Connectivity. 

I. Introduction 

Wireless sensor networks (WSNs) are distributed collections 
of sensors with limited capabilities for computations and 
wireless communications. It is envisioned IT] that WSNs 
will be used in a wide range of applications areas such as 
healthcare (e.g. patient monitoring), military operations (e.g., 
battlefield surveillance) and homes (e.g., home automation and 
monitoring). These WSNs will often be deployed in hostile 
environments where communications can be monitored, and 
nodes are subject to capture and surreptitious use by an adver- 
sary. Under such circumstances, cryptographic protection will 
be needed to ensure secure communications, and to support 
functions such as sensor-capture detection, key revocation and 
sensor disabling. 

Unfortunately, many security schemes developed for general 
network environments do not take into account the unique 
features of WSNs: Public key cryptography is not feasible 
computationally because of the severe limitations imposed on 
the physical memory and power consumption of the individual 
sensors. Traditional key exchange and distribution protocols 
are based on trusting third parties, and this makes them inad- 
equate for large-scale WSNs whose topologies are unknown 
prior to deployment. We refer the reader to the papers [0, 



ifTTI . II20I for discussions of the security challenges in WSN 
settings. 

Random key predistribution schemes were introduced to 
address some of these difficulties. The idea of randomly 
assigning secure keys to sensor nodes prior to network de- 
ployment was first introduced by Eschenauer and Gligor lUTI . 
Since then, many competing alternatives to the Eschenauer 
and Gligor (EG) scheme have been proposed; see for a 
detailed survey of various key distribution schemes for WSNs. 
With so many schemes available, a basic question arises as to 
how they compare with each other Answering this question 
passes through a good understanding of the properties and 
performance of the schemes under consideration, and this can 
be achieved in a number of ways. The approach we use here 
considers random graph models naturally induced by a given 
scheme, and then develops the scaling laws corresponding to 
desirable network properties, e.g., absence of secure nodes 
which are isolated, secure connectivity, etc. This is done 
with the aim of deriving guidelines to dimension the scheme, 
namely adjust its parameters so that these properties occur 
with high probability as the number of nodes becomes large. 

To date, most of the efforts along these lines have been 
carried out under the assumption of full visibility according 
to which sensor nodes are all within communication range 
of each other; more on this later: Under this assumption, the 
EG scheme gives rise to a class of random graphs known 
as random key graphs; relevant results are available in the 
references |I3], JS), ifTTI . ifTSl . Il24l . The q-composite scheme 
Q, a simple variation of the EG scheme, was investigated 
by Bloznelis et al. ||4| through an appropriate extension of 
the random key graph model. Recently, Yagan and Makowski 
have analyzed various random graphs induced by the random 
pairwise key predistribution scheme of Chan et al. Q; see the 
conference papers ||251 . ll26l . 

To be sure, the full visibility assumption does away with 
the wireless nature of the communication medium supporting 
WSNs. In return, this simplification makes it possible to focus 
on how randomization of the key distribution mechanism alone 
affects the establishment of a secure network in the best of 
circumstances, i..e., when there are no link failures. A common 
criticism of this line of work is that by disregarding the 
unreliability of the wireless links, the resulting dimensioning 
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guidelines are likely to be too optimistic: In practice nodes 
will have fewer neighbors since some of the communication 
links may be impaired. As a result, the desired connectivity 
properties may not be achieved if dimensioning is done 
according to results derived under full visibility. 

In this paper, in an attempt to go beyond full visibility, 
we revisit the pairwise key predistribution scheme of Chan et 
al. Q under more realistic assumptions that account for the 
possibility that communication links between nodes may not 
be available - This could occur due to the presence of physical 
barriers between nodes or because of harsh environmental 
conditions severely impairing transmission. To study such 
situations, we introduce a simple communication model where 
channels are mutually independent, and are either on or off. An 
overall system model is then constructed by intersecting the 
random graph model of the pairwise key distribution scheme 
(under full visibility), with an Erdos-Renyi (ER) graph model 
Is). For this new random graph structure, we establish zero- 
one laws for two basic (and related) graph properties, namely 
graph connectivity and the absence of isolated nodes, as the 
model parameters are scaled with the number of users - We 
identify the critical thresholds and show that they coincide. To 
the best of our knowledge, these full zero-one laws constitute 
the first complete analysis of a key distribution scheme under 
non-full visibility - Contrast this with the partial results by Yi 
et al. 1281 for the absence of isolated nodes (under additional 
conditions) when the communication model is the disk model. 

Although the communication model considered here may 
be deemed simplistic, it does permit a complete analysis 
of the issues of interest, with the results already yielding a 
number of interesting observations: The obtained zero-one 
laws differ significantly from the corresponding results in 
the full visibility case IZSl . Thus, the communication model 
may have a significant impact on the dimensioning of the 
pairwise distribution algorithm, and this points to the need 
of possibly reevaluating guidelines developed under the full 
visibility assumption. Furthermore, simulations suggest that 
the zero-one laws obtained here for the on/off channel model 
may still be useful in dimensioning the pairwise scheme under 
the popular, and more realistic, disk model [[TZj . 

We also compare the results established here with well- 
known zero-one laws for ER graphs [5|. In particular, we 
show that the connectivity behavior of the model studied here 
does not in general resemble that of the ER graphs. The 
picture is somewhat more subtle for the results also imply 
that if the channel is very poor, the model studied here indeed 
behaves like an ER graph as far as connectivity is concerned. 
The comparison with ER graphs is particularly relevant to 
the analysis of key distribution schemes for WSNs: Indeed, 
connectivity results for ER graphs have often been used in 
the dimensioning and evaluation of key distribution schemes, 
e.g., see the papers by Eschenauer and Gligor llTTI . Chan 
et al. |7l and Hwang and Kim flSl . There it is a common 
practice to assume that the random graph induced by the 
particular key distribution scheme behaves like an ER graph 
(although it is not strictly speaking an ER graph). As pointed 



out by Di Pietro et al. fS] such an assumption is made without 
any formal justification, and subsequent efforts to confirm its 
validity have remained limited to this date: The EG scheme 
has been analyzed by a number of authors ill], |I8], lITSl . 11241 . 
and as a result of these efforts it is now known that the ER 
assumption does yield the correct results for both the absence 
of isolated nodes and connectivity under the assumption of full 
visibility. On the other hand the recent paper |l25l shows that 
the ER assumption is not valid for the pairwise key distribution 
of Chan et al. Q; see Section [V-AI for details. 

The rest of the paper is organized as follows: In Section |ll] 
we give precise definitions and implementation details of the 
pairwise scheme of Chan et al. while Section |lll] is devoted 
to describing the model of interest. The main results of the 
paper, namely Theorem 14.11 and Theorem 14.21 are presented 
in Section |IV] with an extensive discussion given in Section 
W\ The remaining sections, namely Sections |VI] through IXIIII 
are devoted to establishing the main results of the paper 

A word on notation and conventions in use: All Umiting 
statements, including asymptotic equivalences, are understood 
with n going to infinity. The random variables (rvs) under 
consideration are all defined on the same probability triple 
{il,J-',F). Probabilistic statements are made with respect to 
this probability measure P, and we denote the corresponding 
expectation operator by E. Also, we use the notation =st to 
indicate distributional equality. The indicator function of an 
event E is denoted by 1 [E]. For any discrete set S we write 
5*1 for its cardinality. Also, for any pair of events E and F 
we have 

1[E U F] = 1[E] + 1[F] - 1[E n F] . (1) 
II. Implementing pairwise key distribution 

SCHEMES 

Interest in the random pairwise key predistribution scheme 
of Chan et al. jT] stems from the following advantages over 
the EG scheme: (i) Even if some nodes are captured, the 
secrecy of the remaining nodes is perfectly preserved; (ii) 
Unlike earlier schemes, this pairwise scheme enables both 
node-to-node authentication and quorum-based revocation. 

As in the conference papers ||25l , ||26l , we parametrize the 
pairwise key distribution scheme by two positive integers n 
and K such that K < n. There are n nodes, labelled i = 
1, . . . ,n, with unique ids Idi, . . . , Id„. Write TV := {1, . . . n} 
and set Af-i := N — {i} for each i = 1, . . . , n. With node i 
we associate a subset r„ ^ of nodes selected at random from 
Af-i - We say that each of the nodes in Tn.i is paired to node 
i. Thus, for any subset A C A/"-;, we require 

r (V)"' if|^| = ^r 

[ otherwise. 

The selection of r„ ^ is done uniformly amongst all subsets of 
J\f-i which are of size exactly K. The rvs r„.i, . . . , r„ „ are 
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assumed to be mutually independent so that 

n 

P [r„,, ^A,,t^l,...,n]=Y[F [r„,, = A,] 

1=1 

for arbitrary Ai, . . . , An subsets of A/li, . . . ,7V-„, respec- 
tively. 

Once this offline random pairing has been created, we 
construct the key rings J^n,i, ■ ■ ■ ,^n,n, one for each node, 
as follows: Assumed available is a collection of nK distinct 
cryptographic keys {u}i\£, i — l,...,n; £ = 1,...,A'}. 
Fix i = l,...,n and let in^i : Tn.i {1,...,K} denote 
a labeling of r„ For each node j in r„ ^ paired to i, 
the cryptographic key uj^g^^ .(^j-f is associated with j. For 
instance, if the random set r„ i is realized as {ji, . . . , Jk} 
with 1 < ji < ... < jK < n, then an obvious labeling 
consists in in,i{jk) = k for each k = 1,...,K with key 
w^i^. associated with node jk- Of course other labeling are 
possible, e.g., according to decreasing labels or according 
to a random permutation. Finally, the pairwise key w* ^ = 
[Id,;|Idj|a;j|£^ .(j)] is constructed and inserted in the memory 
modules of both nodes i and j. The key cj* is assigned 
exclusively to the pair of nodes i and j, hence the terminology 
pairwise distribution scheme. The key ring of node i is 
the set 



(2) 



If two nodes, say i and j, are within communication range 
of each other, then they can establish a secure link if at least 
one of the events i £ r„ j or j G Tn.j is taking place. Both 
events can take place, in which case the memory modules of 
node i and j both contain the distinct keys a;* and w* ^j. 
Finally, it is plain by construction that this scheme supports 
node-to-node authentication. 

III. The model 

Under full visibility, this pairwise distribution scheme natu- 
rally gives rise to the following class of random graphs: With 
n ~ 2,3,... and positive integer K < n, we say that the 
distinct nodes i and j are K-adjacent, written i j, if and 
only if they have at least one key in common in their key 
rings, namely 



j iff ^nd n I]„j- ^ 



(3) 



Let H(n; K) denote the undirected random graph on the 
vertex set {!,..., n} induced by the adjacency notion (|3]i; 
this corresponds to modelhng the pairwise distribution scheme 
under full visibility. We have 

P[t^Kj] = KiK) (4) 

where A„(A') is the link assignment probability in M{n;K) 
given by 

/ ft- \ 2 
A„(A-) = 1 - ( 1 



2K 



n-1 
K 



1 



(5) 



As mentioned earUer, in this paper we seek to account 
for the possibility that communication links between nodes 
may not be available. To study such situations, we assume a 
communication model that consists of independent channels 
each of which can be either on or off. Thus, with p in (0, 1), 
let {Bij{p),l < i < j < n} denote i.i.d. {0,l}-valued rvs 
with success probability p. The channel between nodes i and 
j is available (resp. up) with probability p and unavailable 
(resp. down) with the complementary probability 1 — p. 

Distinct nodes i and j are said to be B-adjacent, written 
i j, if Bij{p) = 1. The notion of B-adjacency defines 
the standard ER graph G{n;p) on the vertex set {!,... ,n}. 
Obviously, 

The random graph model studied here is obtained by 
intersecting the random pairwise graph H(n; K) with the ER 
graph G{n;p). More precisely, the distinct nodes i and j are 
said to be adjacent, written i ^ j, if and only they are both 
K-adjacent and B-adjacent, namely 



j iff and 

B,,{p) = 1. 



(6) 



The resulting undirected random graph defined on the vertex 
set {!,..., n} through this notion of adjacency is denoted 
HnG(n; A». 

Throughout the collections of rvs {r„ i, . . . , r„ „} and 
{Bij{p), 1 < i < j < n} are assumed to be independent, in 
which case the edge occurrence probability in H n G{n; K,p) 
is given by 



IV. The results 



(7) 



To fix the terminology, we refer to any mapping A' : No ^ 
No as a scaling (for random pairwise graphs) provided it 
satisfies the natural conditions 



Kn < n, 71 = 1, 2, , 



(8) 



Nn 



(0,1) defines a scaling for 



Similarly, any mapping p 
ER graphs. 

To lighten the notation we often group the parameters K 
and p into the ordered pair 6 = {K,p). Hence, a mapping 
6* : No — > No X (0, 1) defines a scaling for the intersection 
graph H n Gr{n; 9) provided the condition (O holds on the 
first component. 

The results will be expressed in terms of the threshold 
function t : [0, 1] — [0, 1] defined by 



1 







i°g(i-p) 



if p = 
if < p < 1 



if p = 1. 



(9) 



It is easy to check that this threshold function is continuous 
on its entire domain of definition; see Figure |3] 
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A. Absence of isolated nodes 

The first result gives a zero-one law for the absence of 
isolated nodes. 

Theorem 4.1: Consider scalings K : Nq Nq and p : 
No (0, 1) such that 

Pn (2Kn - ^ clogn, n = 1, 2, . . . (10) 



n — 1 ^ 

for some c > 0. 7f lim„^txj Pn 
we have 



p* for some p* in [0,1], then 



lim I 

n— >-oo 



H n G(n; 6'„) contains 
no isolated nodes 

if c< t{p*) 

1 ifc>T{p*). 



(11) 



The condition {TO^ on the scaling No (0,1) x No will 
often be used in the equivalent form 



Pn 2K„ 



n - 1 

with the sequence c : No 



CrJogn, n = l,2, 



(12) 



satisfying lim 



B. Connectivity 

An analog of Theorem 14.11 also holds for the property of 
graph connectivity. 

Theorem 4.2: Consider scalings K ; No ^ No and p : 
No (0, 1) such that O holds for some c > 0. If 
lim„_j.ooPn = P* for some p* in [0, 1], then we have 



lim P [H n 6'„) is connected] 

n— >oo 

if c< t{p*) 

1 if c > t{p*) 



(13) 



where the threshold t{p*) is given by @. 

Comparing Theorem 14.21 with Theorem 14.11 we see that 
the class of random graphs studied here provides one more 
instance where the zero-one laws for absence of isolated 
nodes and connectivity coincide, viz. ER graphs ||5|, random 
geometric graphs lfT9l or random key graphs |l3], [ITSI . [E4l . 

A case of particular interest arises when p* > Q since 
requiring (fTOl i now amounts to 



2if„ 



-^logn 

p* 



(14) 



for some c > 0. Any scaling X : No — ?• Nq which behaves 
like (O must necessarily satisfy Kn = o(n), and it is easy 
to see that requiring (fTOl i is equivalent to 



t log n 



(15) 



for some t > with c and t related by t = si*"- With 
this reparametrization. Theorem 14. 1 1 and Theorem 14.2 1 can be 
summarized in the following simpler form: 



Theorem 4.3: Consider scalings K : No ^ No and p : 
No (0,1) such that lim„_>ooPn = p* > 0. Under the 
condition ifJJb for some t > 0, we have 

lim P [H n G(n; d„ ) contains no isolated nodes] 
= lim P [H n G(?i; 6l„) is connected] 
if t< ?{p*) 



1 if t> ?{p*) 



where we have set 

t{p) 



t{p) 



1 



2p p- log(l - p) ' 



(16) 



0<p<l. (17) 



This alternate formulation is particularly relevant for the 
case Pn = p* (in (0, 1)) for all n — 1,2,..., which captures 
situations when channel conditions are not affected by the 
number of users. Such simplifications do not occur in the more 
realistic case p* = which corresponds to the situation where 
channel conditions are indeed influenced by the number of 
users in the system - The more users in the network, the more 
likely they will experience interferences from other users. 

We now present numerical results that verify ([16). In all 
the simulations, we fix the number of nodes at rt = 200. 
We consider the channel parameters p ~ 0.2, p = 0.4, 
p = 0.6, p = 0.8, and p = 1 (the full visibility case), while 
varying the parameter K from 1 to 25. For each parameter pair 
{K,p), we generate 500 independent samples of the graph 
EI n G{n;K,p) and count the number of times (out of a 
possible 500) that the obtained graphs i) have no isolated 
nodes and ii) are connected. Dividing the counts by 500, we 
obtain the (empirical) probabilities for the events of interest. 
The results for connectivity are depicted in Figure [T] where 
the curve fitting tool of MATLAB is used. It is easy to check 
that for each value of p ^ 1, the connectivity threshold 
matches the prescription (fTSI l. namely K = t{p) log n. It is 
also seen that, if the channel is poor, i.e., if p is close to 
zero, then the required value for K to ensure connectivity 
can be much larger than the one in the full visibility case 
p = 1. The results regarding the absence of node isolation are 
depicted in Figure |2] For each value of p 7^ 1, Figure |2] is 
indistinguishable from Figure [T] with the difference between 
the estimated probabilities of graph connectivity and absence 
of isolated nodes being quite small, in agreement with (fTSI l. 

V. Discussion and comments 
A. Comparing with the full-visibility case 

At this point the reader may wonder as to what form would 
Theorem l4.2l take in the context of full visibility- In the setting 
developed here this corresponds to p = 1 so that EI n G(n; 9) 
coincides with EI(n; K)\ see the curve for p = 1 in Figure [U. 
Relevant results for this case were obtained recently by the 
authors in 1251 . 
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Fig. 1. Probability that HnG(n; K, 'p) is connected as a function of K for 
p = 0.2, p = 0.4, p = 0.6, p = 0.8 and p = 1 with n = 200. 




Fig. 2. Probability that IHInG(n; K,p) has no isolated nodes as a function 
of is: for p = 0.2, p = 0.4, p = 0.6, p = 0.8 and p = 1 with n = 200. 
This figure clearly resembles Figure [T]for all p 1. 



Theorem 5.1: For any K a positive integer, it holds that 

if K =1 



The case where the parameter K is scaled with n is an easy 
corollary of Theorem 15.11 

Corollary 5.2: For any scaling K : ^ No such that 
Kn > 2 for all n sufficiently large, we have the one-law 

lim P[IHI(n;X„) is connecte(l\ = 1. 



Each node in H(n; K) has degree at least K, so that no node 
is ever isolated in H(ri; K). This is in sharp contrast with the 
model studied here, as reflected by the full zero-one law for 
node isolation given in Theorem 14.11 

Theorem 15.11 and its Corollary 15.21 together show that very 
small values of K suffice to ensure asymptotically almost sure 
(a.a.s.) connectivity of the random graph H(n; K). However, 
these two results cannot be recovered from Theorem 14. 2 I whose 
zero-one laws are derived under the assumption p„ < 1 for all 
n = 1, 2, . . .. Furthermore, even if the scaling p : No (0, 1) 
were to satisfy lim„^oo Pn = 1, only the one-laws in Theorem 
14.31 remain since t{p*) = (and t{p*) — 0) at — 1. 
Although this might perhaps be expected given the aforemen- 
tioned absence of isolated nodes in W{n; K), the one-laws for 
both the absence of isolated nodes and graph connectivity in 
H n 9) still require conditions on the behavior of the 
scaling : No — > No, namely ( fTSl l (whereas Corollary 15.21 
does not). 

B. Comparing M n 9) with ER graphs 

In the original paper of Chan et al. Q (as in the reference 
ifOl). the connectivity analysis of the pairwise scheme was 
based on ER graphs fSl - It was assumed that the random 
graph induced by the pairwise scheme under a communication 
model (taken mostly to be the disk model lfT2l ) behaves like an 
ER graph; similar assumptions have been made in ifTTI . llT3l 
when discussing the connectivity of the EG scheme. However, 
this assumption was made without any formal justification. 
Recently we have shown that the full visibility model EI(n; K) 
has major differences with an ER graph. For instance, the 
edge assignments are (negatively) correlated in H(n; K) while 
independent in ER graphs; see l25l for a detailed discussion 
on the differences of H(7i; A') and G{n;p). It is easy to verify 
that the edge assignments in H n G(ri; 6) are also negatively 
correlated; see Section |IX] Therefore, the models EI(?t,; K) 
and H fl 9) cannot be equated with an ER graph, and 
the results obtained in ||251 and in this paper are not mere 
consequences of classical results for ER graphs. 

However, /orma/ similarities do exist between EI n G(n; 9) 
and ER graphs. Recall the following well-known zero-one law 
for ER graphs; For any scaling p : Nq ^ [0, 1] satisfying 

log n 
Pn ^ c 

n 

for some c > 0, it holds that 



lim P [M{n; K) is connected^ = 



1 ifK > 2. 



lim V [G{n;p„) is connected] = 
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P 

Fig. 3. t{p) vs p. Clearly t(p*) = 1 only if limrn-ooPn = = 0. 

On the other hand, the condition (fTOl i can be rephrased more 
compactly as 

log n 

p„A„(X„)~c , c>0 

n 

with the results ( fTTT i and ( fT3] l unchanged. Hence, in both ER 
graphs and H n 9), the zero-one laws can be expressed 
as a comparison of the probability of link assignment against 
the critical scaling iHIii- this is also the case for random 
geometric graphs |[T9l . and random key graphs ||3], llTSl . Il24l . 
But the condition c > t(p*) that ensures a.a.s. connectivity in 
EI n 6) is not the same as the condition c > 1 for a.a.s. 
connectivity in ER graphs; see Figure|3] Thus, the connectivity 
behavior of the model H n G(n; 6) is in general different 
from that in an ER graph, and a "transfer" of the connectivity 
results from ER graphs cannot be taken for granted. Yet, the 
comparison becomes intricate when the channel is poor: The 
connectivity behaviors of the two models do match in the 
practically relevant case (for WSNs) lim„^ooPn = since 
r(0) = 1. 

C. A more realistic communication model 

One possible extension of the work presented here would 
be to consider a more reahstic communication model; e.g., 
the popular disk model lfT2l which takes into account the geo- 
graphical positions of the sensor nodes. For instance, assume 
that the nodes are distributed over a bounded region T) of the 
plane. According to the disk model, nodes i and j located at 
Xi and Xj, respectively, in V are able to communicate if 

II Xi - Xj II < p (18) 

where p > is called the transmission range. When the node 
locations are independently and randomly distributed over the 




Fig. 4. Probability that H n G(n; K, p) is connected as a function of K. 
The number of nodes is set to n = 200 and p is given by irp^ = p. 

region V, the graph induced under the condition ([18) is known 
as a random geometric graph llT9l , thereafter denoted G{n; p). 

Under the disk model, studying the pairwise scheme of 
Chan et al. amounts to analyzing the intersection of H(n; K) 
and G(n;p), say H n G(n; if, p). A direct analysis of this 
model seems to be very challenging; see below for more 
on this. However, limited simulations already suggest that 
the zero-one laws obtained here for Hn(G(n;A', p) have an 
analog for the model H n K, p). To verify this, consider 
200 nodes distributed uniformly and independently over a 
folded unit square [0,1]^ with toroidal (continuous) boundary 
conditions. Since there are no border effects, it is easy to check 
that 

¥[\\xi-Xj\\< p]^TTp^, i^j, i,j = l,2,...,n. 

whenever p < 0.5. We match the two communication models 
p) and (G(n; p) by requiring Tip^ = p. Then, we use the 
same procedure that produced Figure [T]to obtain the empirical 
probability that H n p) is connected for various values 

of K and p. The results are depicted in Figure |4] whose 
resemblance with Figure [1] suggests that the connectivity 
behaviors of the models H n G(7i; K,p) and H n K, p) 
are quite similar. This raises the possibility that the results 
obtained here for the on/off communication model can also 
be used for dimensioning the pairwise scheme under the disk 
model. 

A complete analysis of H n K, p) is likely to be very 
challenging given the difficulties already encountered in the 
analysis of similar problems. For example, the intersection of 
random geometric graphs with ER graphs was considered in 
II2I, ll28l . Although zero-one laws for graph connectivity are 



Fig. 5. An instantiation of ER grapli G(n;p) with n = 50 and p = O.2.- 
The grapli is connected. 



Fig. 6. An instantiation of EI(n; K) with n = 50 and K = 5.- The graph 
is connected. 



available for each component random graph, the results for the 
intersection model in 121, ll28l were limited only to the absence 
of isolated nodes; the connectivity problem is still open for that 
model. Yi et al. |28ll also consider the intersection of random 
key graphs with random geometric graphs, but these results 
are again limited to the property of node isolation. To the best 
of our knowledge. Theorem 14.21 reported here constitutes the 
only zero-one law for graph connectivity in a model formed 
by intersecting multiple random graphs! (Except of course the 
trivial case where an ER graph intersects another ER graph.) 

D. Intersection of random graphs 

When using random graph models to study networks, sit- 
uations arise where the notion of adjacency between nodes 
reflects multiple constraints. This can be so even when dealing 
with networks other than WSNs. As was the case here, such 
circumstances call for studying models which are constructed 
by taking the intersection of multiple random graphs. However, 
as pointed out earlier, the availability of results for each 
component model does not necessarily imply the availability 
of results for the intersection of these models; see the examples 
provided in the previous section. 

Figures |5]|7] can help better understand the relevant issues 
as to why this is so: Figure |5] provides a sample of an ER 
graph G{n,p) with n = 200 and p = 0.2. As would be 
expected from the classical results, the obtained graph is very 
densely connected. Similarly, Figure |6] provides a sample of 
the pairwise random graph Ell(n; A') with n = 200 and K = 5. 
In line with Theorem l5.ll the obtained graph is connected. On 
the other hand, the graph formed by intersecting these graphs 
turn out to be disconnected as shown in Figure [T] 




Fig. 7. The intersection H n G(n; 9) of the graphs in Figure [s] and Figure 
|6]- The graph is disconnected as the marked nodes form a component! 



To drive this point further, consider the constant parameter 
case for the models M{n;K) and G{n;p), a case which 
cannot be recovered from either Theorem 14.1 l or Theorem I4.2J 
Nevertheless, Theorem 15.11 yields 

lim P [H(n; K) is connected] = 1, K >2 

n—^oo 



g 



while it well known fS) that 

lim P [G{n;p) is connected] = 1. < p < 1. 
However, it can be shown that 

lim P [H n 9) contains no isolated nodes] = (19) 



whence 



lim ] 

n— f oo 



" n G{n; 9) is connected] = 



(20) 



for the same ranges of values for p and K; for details see 
the discussion at the end of Section 1X1 This clearly provides 
a non-trivial example (one that is not for an ER intersecting 
an ER graph) where the intersection of two random graphs is 
indeed a.a.s. not connected although each of the components 
is a.a.s. connected. 

VI. A PROOF OF TheoremI4.1I 



We prove Theorem 14.11 by the method of first and second 
moments lfT4l p. 55] applied to the total number of isolated 
nodes in H n G{n; 9). First some notation: Fix n = 2,3,... 
and consider 9 = {K,p) with p in (0, 1) and positive integer 
K such that K < n. With 

XnA^) ■= 1 [Node i is isolated in H n G{n\ 9)] 

for each i = the number of isolated nodes in 

H n G(n; 9) is simply given by 



Ini9) ■.= J2xnA(^)- 



The random graph H fl G{n; 9) has no isolated nodes if and 
only if /„(6') ^ 0. 

The method of first moment llT4l Eqn (3.10), p. 55] relies 
on the well-known bound 



l-E[/„(0)] <P[/„(0) =0] 



(21) 



while the method of second moment lfT4l Remark 3.1, p. 55] 
has its starting point in the inequality 

The rvs Xn,i{&)i ■ ■ ■ ,Xn,n{(^) being exchangeable, we find 

E[/„(^^)] =nE[x„,iW] (23) 

and 

E [U9f] = nE [xnjrn + n{n - 1)E [xn,i(^)x„.2 W] 
by the binary nature of the rvs involved. It then follows that 
E [U9)^] _ 1 



E[In{9)f nE[xn.im 

E[x„4(g)Xn,2W] 
n ■ (E[x„,i(^)])' 



(24) 



From and (|23l l we see that the one-law 

linin^co P [In{On) = 0] = 1 will be established if we 
show that 

lim nE[x„,i(^n)] =0. (25) 

It is also plain from (|22] | and ( l24l i that the zero-law 

lim„^oo P [In{On) = 0] = holds if 



lim nE [xn, 1(6*70] = oo 



(26) 



E[x„,(^„)x..y]Ul. (27) 

r.^oo'^l (E[x„,l(^n)]) /" 



and 



The proof of Theorem 14. 1 [ passes through the next two tech- 
nical propositions which establish (IZSl l. ( |26] | and ( |27] | under 
the appropriate conditions on the scaling 6* : Nq — > Nq x (0, 1). 

Proposition 6.1: Consider scalings A' : No — > No and p : 
Nq (0, 1) such that dTol ) holds for some c > 0. Assume also 
that lim„_>oo Pn ~ P* exists. Then, we have 

( if c> Tip*) 

\mY^nE[xn,i{Gn)] = \ (28) 
" °° [ oo ifc < t{p*) 

where the threshold t{p*) is given by @. 

A proof of Proposition 16.11 is given in Section IVIIII 

Proposition 6.2: Consider scalings if : No — >■ No and p : 
Nq — > (0, 1) such that dTOl ) holds for some c > 0. Assume also 
that lim„_^oo Pn = P* exists. Then, we have ^27\) whenever 
p* < 1. 

A proof of Proposition 16.21 can be found in Section 
1X1 To complete the proof of Theorem 14.1] pick a scaling 
9 : No -> No X (0, 1) such that ([TOl) holds for some 
c > and \im.n^ooPn = P* exists. Under the condition 
c > t(p*) we get dZST l from Proposition 16.11 and the one- 



law lim 



Ir, 



0] = 1 follows. Next, assume that 
c < t{p*) - This case is possible only if < 1 since 
t(1) = as seen at When p* < 1, we obtain and (|27l ) 
with the help of Propositions 16.1! and 16.21 respectively. The 
conclusion lim„^oo P [^n(^n) = 0] = is now immediate. 

VII. A PREPARATORY RESULT 

Fix ?i = 2, 3, . . . and consider 9 = {K,p) with p in (0, 1) 
and positive integer K such that K < n. Under the enforced 
assumptions, for alH = 1, . . . , n, we easily see that 



E[x,m(0)]-E[(1-p)^ 



(29) 



where Dn,i denotes the degree of node i in M{n;K). Note 
that 

71 

Am = k+ i[*er„,,]. (30) 

j=ij^r„,,u{i} 

By independence, since 

|{j = 1, . . . , n : j ^ r„,, U {i}}\ ^n-K-l, 
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the second term in ( [3Q] l is a binomial rv with n — K —1 trials 
and success probability given by 



K 



whence 



(31) 



(32) 



The proof of Proposition 16. 1 l uses a somewhat simpler form 
of the expression (|32] i which we develop next. 

Lemma 7.1: Consider scalings A' : No — !• No andp : Nq — s- 
(0,1) such that (W^ holds for some c> 0. It holds that 



n=l,2,... (33) 



with 



an ■■= (1 - C„) logn + KniPn + log(l - Pn)) (34) 

where the sequence c : Nq — > K is the one appearing in the 
form {121) of the condition (fJOl ). 

In what follows we make use of the decomposition 



log(l - x) = -X - -^(x), < a; < 1 (35) 



with 



1-i 



-dt 



on that range. Note that 



lim 



^{x) 1 



where the last step used the form (fT2] i of the condition JTOl i 
on the scaling. Reporting this calculation into the expression 
for /?„ we find 



fin = an-{n- Kn - 1)* 



PnKn 



Lemma PTTI will be established if we show that 
lim {n - A-„ - 1)^ (P^) = 0. 
To that end, for each n = 2, 3, ... we note that 



(37) 



PnKn < Pn 2K, 



< 2pnKn 



n — 1 

since A'„ < n. The condition ( fT2l i implies 

Y logn < PnKn < Cn logn, 
and it is now plain that 



(38) 



lim = and lim (n - A'„ - 1) - 0. 



n — 1 



(n - 1)2 



Invoking the behavior of 'i'{x) at a: = mentioned earlier, we 
conclude from these facts that 



lini_ ( [n ~ Kn - I)— —2 



(n - 1)2 



n~l 




= 0. (39) 



This establishes ( |37] i and the proof of Lemma 17.11 is 
completed. ■ 



Proof. Consider a scaling 6* ; No ^> No x (0, 1) such that 
dTol i holds for some c > and assume the existence of the 
limit lim„_>ooP?i = P*- Replacing 6 by On in (|32] | for each 
ri = 2, 3, . . . we get 



with [3n given by 

/3„ = logn + Kn log(l - p„) - 7„ 

with 



(36) 



7n 



-(n-/^„-l)log 1- 



PnKn 

n - 1 



The decomposition ( |35] l now yields 

7„ := (n-A-„-l)('^ + *^^"^^" 



n — 1 



n - 1 



= ( 1 - KnPn + {n- Kn 1)^ 



n-1 



n — 1 



-KnPn + 2 - 



n — 1 



- /Y„ - 1)* 



KnPn 
PnKn 

n- 1 



-KnPn + Cn logU + {u - Kn - 1)^* 



PnKn 

n ~ 1 



VIII. A PROOF OF PropositionI6.1 



In view of Lemma 17. 1[ Proposition 16. II will be established 
if we show 



lim ar, 

n— >oo 



-oo if c > r(j3*) 
-oo if c < t{p*). 



(40) 



To see this, first note from (135b that for each n = 1,2, 
we have p„ + log(l — Pn) < and the lower bound in 
implies 



an < (1 -c„)logri + c„ 



log n 
2p,i 



(p„ +l0g(l -pn)) 



Letting n go to infinity in this last expression, we get 

lim„^oo ctn = —oo whenever 

2 



> lim 



since lim„_>.oo c„ = c. 



1 - 



log(l-p„) 



(42) 



Next, we show that if c < t(j)*), then lim„^oo Qfn = +oo. 
We only need to consider the case < < 1 since r(l) = 
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and the constraint c < t(1) is vacuous. We begin by assuming 
p* = 0, in which case for each n = 2, 3, . . ., we have 

Ctn = (1 - C„) logn + Kn{pn + {-Pn ^ *(PrO)) 

= (1 - c„) log n - ( — — ■ KnPn 



Pn 



> (1 - c„) log n - Cn log n • [ — ^j— ) p„ 

Pn 



log n • 1 - c„ 1 



pI 



Pr 



(43) 



with the inequality following from the upper bound in (l38T l. 
Let n grow large in the last expression. Since we have assumed 
lim„_^oo Pn = 0, we get 



lim Pn 



■^{Pn) 



Pn 



0, 



and the desired conclusion lim„_j.oo Un 



~oo is obtained 



whenever c < 1 = t(0) upon using lim„^oo c„ = c. 

Finally we assume < < 1. For each e > 0, there exists 
a finite positive integer n*{e) such that p„ > (1 — e)p* when 
n > n*{e). On that range the upper bound in (l3Ft yields 



A' < 



loeri, 



(1 - e)p* 

whence the conclusions ~ o{n) and 



Pn 2K 



71-1 



2KnPn + 0(1) 



follow. Comparing this last fact against the lefthand side of 
(fni yields 



so that 



KnPn = Y logn + 0(1), 



ATnPn Y logn. 



(44) 



From (|34] l it follows that 

"- = (l-c„) + 



logn 



Pn J log n 



for all n sufficiently large. Letting n go to infinity in this last 
expression and using (|44] | with the earlier remarks, we readily 
conclude 

lin, -^ = (l-c) + ^('l + Ml^Ul__£_ 
n-)-oo logn 2 \ P* J tkP*) 

where the last step follows by direct inspection. It is now 
clear that lini„^oo = oo when c < t{p*) with < < 1. 
The proof of Proposition 16.11 is now completed. ■ 



IX. Negative dependence and consequences 

Fix positive integers n ~ 2,3,... and K with K < n. 
Several properties of the {0, l}-valued rvs 



1 [j e r„,,] 



i 7^ j 
i, j = 1, . . . ,n 



(45) 



and 



1 [j e r„,, V I e r„ j-] , , , ' ^ ^ (46) 



i 7^ j 
i,j = l,...,n 



will play a key role in some of the forthcoming arguments. 

A. Negative association 

The properties of interest can be couched in terms of 
negative association, a form of negative correlation introduced 
to Joag-Dev and Proschan ITSl . We first develop the needed 
definitions and properties: Let {X\, A € A} be a collection 
of R-valued rvs indexed by the finite set A. For any non- 
empty subset A of A, we write Xa to denote the R'^'-valued 
Xa = (Xx, X e A). The rvs {Xx, A € A} are then said 
to be negatively associated if for any non-overlapping subsets 
A and i? of A and for any monotone increasing mappings 
(/3 : mI"^I — > K and t/i : M'^' — > R, the covariance inequality 

E [^{XaWXb)] < E MXa)] E mXe)] (47) 

holds whenever the expectations in ( |47] | are well defined and 
finite. Note that Lp and ip need only be monotone increasing 
on the support of Xa and Xg, respectively. 

This definition has some easy consequences to be used 
repeatedly in what follows: The negative association of 
{Xx, A G A} implies the negative association of the collection 
{Xx, A e A'} where A' is any subset of A. It is also well 
known ITSl P2, p. 288] that the negative association of the rvs 
{Xx, A S A} implies the inequality 



E 



n 



LxeA 



< 



l[E[fxiXx 



(48) 



xeA 



where A is a subset of A and the collection {fx, A e A} 
of mappings M — > are all monotone increasing; by non- 
negativity all the expectations exist and finiteness is moot. 

We can apply these ideas to collections of indicator rvs, 
namely for each A in A, Xx = 1 [Ex] for some event Ex- From 
the definitions, it is easy to see that if the rvs {1 [Ex] , A G A} 
are negatively associated, so are the rvs {1 [£'^] , A G A}. 
Moreover, for any subset A of A, we have 



'[Ex, XeA]< 1[f[Ex]. 



(49) 



This follows from (l48T l by taking fx{x) = a;"*" on R for each 
A in A. 
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B. Useful consequences 

A key observation for our purpose is as follows: For each 
i ~ 1, . . . ,n, the rvs 

{i[je r„,,] , J G A/-^ J (50) 

form a collection of negatively associated rvs. This is a 
consequence of the fact that the random set r„ ^ represents 
a random sample (without replacement) of size K from Af-f, 
see lITSl Example 3.2(c)] for details. 

The n collections (|50] | are mutually independent, so that by 
the "closure under products" property of negative association 
ifBl P7, p. 288] GOl p. 35], the rvs ^ also form a collection 
of negatively associated rvs. 

Hence, by taking complements, the rvs 



1 b i r„,,;] 



« 7^ J 

j,J = 1, • • 



(51) 



also form a collection of negatively associated rvs. With 
distinct i, j = 1, . . . , n, we note that 

1 {% i r„j-, J i r„,,] = f{i[ii r„,,] , i [j ^ r„,,]) (52) 



with mapping / 



given by f{x,y) = for 



all X, y in R. This mapping being non-decreasing on 
follows IB] P6, p. 288] that the rvs 

i,j = l,...,n 



(53) 



are also negatively associated. Taking complements one more 
time, we see that the rvs ( |46] | are also negatively associated. 

For each fc = 1, 2 and j ^ 3, . . . ,n, we shall find it useful 
to define 



and 



Under the enforced assumptions, we have hn.ziO) = 

bn,n{0) = K{0) and U„,3,i(6') = ... = lt„,ri,l(^) = 
Un,S,2{6) = . . . = Un,n,2{d) = Un{0)- 

Before computing either one of the quantities Un{9) and 
bn{0), we note that 



bn{0) < Un{0f 



(54) 



This is a straightforward consequence of the negative associ- 
ation of the rvs ( |45] | - In (|47] |. with A and B singletons, use 
the increasing functions iy9,-0:R^M;a-— !>— (1— p)^. 
Using dSTJ we get 



un{e) 



l-p 



K 

n - 1 
K 



K 

n- 1 



n — 1 



(55) 



An expression for bn{0) is available but will not be needed 
due to the availability of ( |54] |. 



X. A PROOF OF Proposition I6.2| 

As expected, the first step in proving Proposition 
consists in evaluating the cross moment appearing in the 
numerator of (l27t . Fix n = 2, 3, . . . and consider 6 = {K,p) 
with p in (0, 1) and positive integer K such that K < n. 
Define the No-valued rvs Bn{9) and C/„(6') by 



:= ^ 1 [j ^ r„,i] 1 [j ^ r„,2] 

3=3 



and 



J2i[3^ r„,i] 1 b 6 r„,2] 



(56) 



(57) 



i=3 



i=3 

Proposition 10.1: Fix ?i = 2, 3, . . .. For any p in (0, 1) and 
positive integer K such that K < n, we have 

E[Xn.l{0)Xn.2m (58) 



(1 -p)'^E 



(1 _p)i[2Gr„,i,ier„,2] 

where the rvs En {0) and Un [8) given by l56l ) and (f57l ). respec- 
tively. 

A proof of Proposition 110.11 is available in Appendix |A] 
Still in the setting of Proposition 110.11 we can use ( l54l i in 
conjunction with (l58T l to get 

nxnAe)xn:2m 



(59) 



It is plain that 

2B„(0) + [/„(0) 



(1 -p)i[2er„,i,ier„,2] 



^ 1 [j ^ r„a] + ^ 1 [.? ^ r„,2 



J=3 



J=3 



= We note that 



E 1 ^ r„,i] 

i=3 



E 1 [.? ^ r„,i] - 1 [2 ^ r„,i] 
(n - 1 - i^) - (1 - 1 [2 e r„,i] 

(n - 2 - X) + 1 [2 G r„,i] 



and 



E 1 [■?■ ^ r„,2] = (n - 2 - X) + 1 [1 e r„,2] 

J=3 

by similar arguments. The expression 

2B„(0) + [/„(0) (60) 
= 2(n - 2 - + 1 [2 € r„a] + 1 [1 G r„,2] 
now follows, and we find 

nXnAe)Xn.2{e)] (61) 
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with 



Rn{0) :=E 



UniO) 



i[2Gr„,i]+i[i6r„ 



(1 _p)i[2er„,i4er„,2: 
Next, with the help of (|32] | and i55[ we conclude that 

E[Xn,l(^)Xn,2(^)] 



< 



2(n-2-A:) 



ii„(0)i[2er„,i]+i[ier„,.]-2 



E 



(1 _p)i[2er„,i,iGr„,2] 

Direct inspection readily yields 

?/„(0)i[2er„,i]+i[i6r„.2]-2 



(1 -p)i[26r„,i4Gr„,2] 



(62) 



(63) 



1 

i-p 



1 - 



1 



pK 



pK 



if 2er„,i,i 6r„,2 



if 2 ^r„,i,i ^r„,2 



otherwise. 



Taking expectation and reporting into (162) we then find 



(E[X„4W]) 
1 



(64) 



< 



_[2er„,i,ier„.2]+ 1-p 

I — p \ n — 1 



K 



I — p \n — 1 



K 



l-p 



K 

n - 1 



(65) 



by a crude bounding argument. 

Now consider a scaling : No ^> No x (0, 1) such that (fTOl i 
holds for some c > and lim„_>.oo Pn = p* < 1. Replace 
6 by On in the bound (l65l l with respect to this scaling. It is 
immediate that ( |27] | will be established if we show that 

2 



1 



and that 



lim 

n->oo 1 — p„ \ n — 1 

n — 1 



lim 1 - p„ 



= 



= 1. 



These limits are an easy consequence of the inequaUties 
by virtue of the fact that lim„^ooPri = p"^ < 1. 



We close with a proof of ( fT9] i: Consider 9 = {K,p) with p 
in (0, 1) and positive integer K. It follows from (|32] | that 



lim E[xn.i{0)] = 



whence lim„^oo E [/n(^)] = oo. It also immediate from 
that 

,^3,pE[x„,lWX.2W]^^^ 
(E[Xn,l(f?)]) 



The arguments outlined in Section |VT] now yield 

lim P[/„(e') = 0] = 0, 

and this establishes (fT9] l. The conclusion ( |20] | immediately 
follows; see discussion at (|66] l. 



XI. A PROOF OF Theorem |4.2| (Part 1) 

Fix n = 2, 3, . . . and consider 6 = {K,p) with p in (0, 1) 
and positive integer K such that K < n. We define the events 



Cn{0) := [HnG(n;6') is connected] 



and 



i{n-e) 



:n(G n; 



contains no isolated nodes] 



If the random graph H n 0) is connected, then it does 
not contain any isolated node, whence Cn{0) is a subset of 
I(n: 6), and the conclusions 



P[Cn{0)]<P[Iin;e)] (66) 



and 



p [Cr^iey] = p [Cniey n i{n; e)] + p [i{n; ey] (6?) 

obtain. 

Taken together with Theorem 14.11 the relations (|66t and 
(|67] | pave the way to proving Theorem 14.21 Indeed, pick a 
scaling 6* : No -> No x (0, 1) such that ([TOll holds for some 
c > and lim„_>oo Pn = P* exists. If c < t(p*), then 
liiTLn^oo P On)] = by the zero-law for the absence 
of isolated nodes, whence lini„_>oo P [Cri(^ri)] ^ with the 
help of dSHl. If c > t(p*), then lini„^oo P [/(n; 6'„)] = 1 
by the one-law for the absence of isolated nodes, and the 
desired conclusion lim„^oo P [C'»i(^?i)] = 1 (or equivalently, 
lim„^oo P [CniOnY] ~ 0) will follow via d67T i if we show the 
following: 

Proposition 11.1: For any scaling : Nq NqX (0, 1) such 



that lim. 
we have 



n— >oo Pn 



p* exists and (fJOl ) hoJds for some c > t{j)*). 



lim P[C„(0„)'n/(7i;0„)] = 0. 



(68) 



The proof of Proposition 111.11 starts below and runs through 
two more sections, namely Sections IXllI and IXlllI The basic 
idea is to find a sufficiently tight upper bound on the proba- 
bility in ( |68] | and then to show that this bound goes to zero 
as n becomes large. This approach is similar to the one used 
for proving the one-law for connectivity in ER graphs p. 
164]. 

We begin by finding the needed upper bound: Fix n = 
2,3,... and consider = {K,p) with p in (0,1) and 
positive integer K such that K < n. For any non-empty 
subset S of nodes, i.e., S C we define the 

graph H n G{n; 0){S) (with vertex set S) as the subgraph of 
H n G(n; 0) restricted to the nodes in S. We also say that S is 
isolated in EI H G(n; 0) if there are no edges (in H H G(n; 0)) 
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between the nodes in S and the nodes in the complement 
5^ = {1, . . . , ri} — 5. This is characterized by 

E„,, n E„j = V B,j{p) =0, ieS, J e S^. 

With each non-empty subset 5* of nodes, we associate 
several events of interest: Let C„(0; S) denote the event that 
the subgraph Mr\G{n;9){S) is itself connected. The event 
Cn{0; S) is completely determined by the rvs {Ki{9), i G 5}. 
We also introduce the event Bn{0; S) to capture the fact that 
S is isolated in H n 9), i.e., 



Bnie;S) 

n E„,, = V B,,{p) =0, i E S, j e S'] . 

Finally, we set 

A^{9;S) ■.^a,{e;S)nB,,{9;S). 

The starting point of the discussion is the following basic 
observation: If H n G(n; 9) is not connected and yet has no 
isolated nodes, then there must exist a subset S of nodes with 
5*1 > 2 such that Mn G{n; 9){S) is connected while S is 
isolated in H n G(n; 9). This is captured by the inclusion 

Cn{9rnl{n;9) C |J (69) 

SCAT: |S|>2 

A moment of reflection should convince the reader that this 
union need only be taken over all subsets 5* of {!,..., n} 
with 2 < \S\ < [^J. A standard union bound argument 
immediately gives 



v[Cn{erni{n; 



SC7V:2<|S|<[-|J 

= E ( E PK(^;5)] |(70) 

where A/",!,,, denotes the collection of all subsets of {1, . . . ,n} 
with exactly r elements. 

For each r = 1 , . . . , n, we simplify the notation by writing 
AnAO) ■■= An{9;{l,...,r}), B„,,(0) := B„(0; {1, . . . , r}) 
and Cn,ri9) := Cn{9; {1, . . . ,r}). With a slight abuse of 
notation, we use Cn{9) for r = n as defined before. Under 
the enforced assumptions, exchangeability yields 

and the expression 

^ ¥[A,,i9;S)] = PK,,.(0)] (71) 

follows since |7V„^r| = (")■ Substituting into ( fTOl i we obtain 
the key bound 

LtJ . X 

P [C„{9r n I{n; 9)]<J2[)P Kx(^)] • (72) 

Consider a scaling : Nq ^> No x (0, 1) as in the statement 
of Proposition ! 1 1.11 Substitute 9 by 0„ by means of this scaling 



in the right hand side of (l72l l. The proof of Proposition 111.11 
will be completed once we show 



LtJ 



lim y 



r=2 



[AnA&n)] = 0. 



(73) 



The means to do so are provided in the next section. 

XII. Bounding probabilities 

Fix n = 2, 3, . . . and consider 9 = {K,p) with p in (0, 1) 
and positive integer K such that K < n. 

A. Bounding the probabilities P [Bn^r{Q)\ 

The following result will be used to efficiently bound the 
probability P[B„,^ (61)]. 

Lemma 12.1: For each r = 2,..., n — 1, we have the 
inequality 



Bn,r(9) 



(74) 



with Un{0) defined by 05l ) and ttie rv ^ given by 

n r 

Kr ■■= E E 1 ^ r„.,] . 

i=r+l 1=1 



(75) 



A proof of Lemma [12. II is available in Appendix iB] The rv 
-E* ^, which appears prominently in (|74] |. has a tail controlled 
through the following result. 

Lemma 12.2: Fix r = 2, . . . , n — 1. For any t in (0, 1) we 
have 



< (1 - t)rX 



n- 1 



(76) 



Proof. Fix n = 2,3,... and consider a positive integer K 
such that K < n. From the facts reported in Section |IX] the 
negative association of the rvs dSOl l implies that of the rvs 
{1 [e e Tn,i] , i = r + 1, . . . , n; ^ = 1, . . . , r}. We ai'e now 
in position to apply the Chemoff-Hoeffding bound to the sum 
drsl ). We use the bound in the form 

P [El, < (1 - t)E [El,]] < e-#«K.^] (77) 

as given for negatively associated rvs in |fTO]. Thm. 1.1, p. 6]. 
The conclusion ( |76] | follows upon noting that 



^[Elr]-= E E^[^er„,,]=r(n-r) 



i=r+l £=1 



n — 1 



as we use 



(EB. 
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B. Bounding the probabilities P [Cn,r{0)\ 

For each r = 2,...,n, let W^Gjr{n\9) stand for the 
subgraph H n G(n; e){S) when 5 = {1, . . . , r}. Also let % 
denote the collection of all spanning trees on the vertex set 

Lemma 12.3: Fix r = 2, . . . ,n. For each T in%., we have 



[rcHnG,(n;0)] < {p\n{K))' 



(78) 



where the notation T C H n indicates that the tree T 

is a subgraph spanning H n Gr {n; 0). 

Since pXn{K) is the probability of link assignment, the 
situation is reminiscent to the one found in ER graphs |j5] 
and random key graphs 1,23 i where in each case the bound 
( ITST i holds with equality. 

Proof. Fix r = 2,3, ... ,n and pick a tree T in %■ Let 
£{T) be the set of edges that appear in T. It is plain that 
T C H n Gr(n, ; 0) occurs if and only if the set of conditions 



and 



holds. Therefore, under the enforced independence assump- 
tions, since = r — 1, we get 

¥[T (imnGr{n;0)] 



Y[ 1 [s„,, n E„,, ^ 

l.j:{l.j}e£{T) 



Y[ 1 e r„j- V J e r„,,] 

_i.j:{l.j}e£(T) 

Y[ P [z e r„ J- V J e r„,,] (79) 

by making use of ( |49l ) with the negatively associated rvs 
( |46] |. The desired result (fTsT l is now immediate from (|5) and 
the relation \£{T)\ ^ r - 1. ■ 



< 



As in ER graphs ||5] and random key graphs 
to the following bound. 

Lemma 12.4: For each r = 2, . . . ,n, we have 



(80) 



Proof. Fix r = 2,...,n. If IHInGr(n;0) is a connected 
graph, then it must contain a spanning tree on the vertex set 
{!,... .r}, and a union bound argument yields 

P [CnAS)] < XI IP ^ ^ G{n; e){S)] . (81) 

By Cayley's formula fT6T| there are r^~'^ trees on r vertices, 
i.e., \Tr\ = r*""^, and dSOl l follows upon making use of (fTSl l. ■ 



XIII. A PROOF OF Proposition II I.IK Part II) 

Consider a scaling : No ^ No x (0, 1) as in the statement 
of Proposition [nH Pick integers i? > 2 and n*{R) > 2{R + 
1) (to be specified in Section lXIII-BI) . On the range n > n*{R) 
we consider the decomposition 

Lf J . ^ 



r=2 ^ ^ r=ii;+l ^ ^ 

and let n go to infinity. The desired convergence (|7Jt will be 
established if we show 



lim 



for eac/j r = 2, 3, . . . and 

LtJ 

lim V 

rt. — ^oo ' ^ 



n— >oo ^ — ' \ r 

r=R+l 



K,r(0„)]=O (82) 



'K,r(^n)]=0. (83) 



We establish (|82] | and (l83T l in turn. Throughout, we make 
use of the standard bounds 



n\ /en 
r / \ r 



r = 1, 



(84) 



for each n = 2, 3, , 

A. Establishing 

Fix r = 2, 3, . . . and consider n = 2, 3, . . . such that r < n. 
Also let = {K,p) with p in (0,1) and positive integer K 
such that K < n. With dTSI ) in mind, for each i = 1, . . . , r, 
we note that 

n 71 r 

^iker^A - X 1 [fc e r„,,] - X 1 [A: G r„,,] 

fe=r+l 



fe=l 



fe=l 



k=l 

we have since |r„ = K. The bounds 

n 

{K-r)+ < J2 1 e < K 
follow, whence 



(85) 



k=r+l 



r{K - r)+ < E*^^ < rK. 
It is also the case that 

r{n — 7' — K)'^ < r{n — r) — i?* 
Reporting these lower bounds into (l74l . we get 



+ 



< (l-p)''^^-'"' •u„(0)'-("-'-^) 



(86) 
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since < p, Un{0) < 1. If we set 
it is now plain that 



= E 



l[C7„,r(^)]P S„,,(0) 



r 



Applying Lemma ri2.4l we find 



(87) 



< 



< 



as we make use of ( |84] i. 
We also note that 



(88) 



(89) 



with 

FAW 



{K — r) log(l — p) — {n — r — K)p 

K 



[K - r) log(l - p) - 1 - 



K 

n — 1 
r - 1 



n — 1 n — 1 



pK 



= [K — r) log(l — p) — p [K — 



n — 1 

X(p + log(l -p)) - rlog(l -p) 



r - 1 
n- 1 



p 2A' - 



n — 1 



n - 1 



(90) 



Now, pick any given positive integer r = 2, 3, . . . and con- 
sider a scaling 6* : No No x (0, 1) such that lim„_>tx> Pn = P* 
exists and (fTOl i holds for some c > t(p*). Replace by 6'„ 
in dSSl ) according to this scaling. In order to establish (l82t it 
suffices to show that 



lim (en)'^ (p„A„(A'„))'' 



(91) 



For n sufficiently large, from (fT2t and (ISFt we first get 



P[A„.,(0)] 



< {eny {pnXn{KAy~^ ■ FnA^A" 

r-1 



loe 



n 



71 — 1 

n 



F, 



" 1 

n — 1 



logn 



K,,(0„)^ (92) 



On the other hand, upon making use of the bounds at 
we find 

FAiOA < ^r„(p„ + log(l-p„))-rlog(l-p„) 



Pn ^K, 



Kl \ r 



PnKn 



n — 1 / n 

= Kn (pn + log(l - pA) - r log(l - pA 
r 

- C„logn+ - PnKn 

n 

< Kn ipn + l0g(l - pA) - Cn log Tl 

T 

-r log(l - pA + - c„ log n 
n 

= p„A„ 1 H Cn logn 

V Pn J 

T 

- r log(l -pn) + - Cn log n 

n 

^ c„ / log(l -pA\ , 

< — log n • IH Cn log n 

2 V Pn J 

T 

- r log(l -pn) + - Cn log n 

n 

Cn A l0g(l-p„)\i 

- rlog(l - Pn) + - c„logn. 

n 

2rp 



= log". -- 



log n 



1 - 



log(l -pA 



Pn 



rpn + - Cn log n 



^ logn 2rpn \ A log(l-p„) 

log n 



H c„logn. 

n 



(93) 



(94) 



As a result, (|90] | implies 
Under the enforced assumptions of Theorem 14.21 we get 

li^ fi_!:L_^Vfi_Mi^ 

n^oo \ 2 \ logn J \ Pn 

= I -A- (^i_ML^ 

2 V 

r(p*) 

and the desired conclusion ( |9T1 i follows upon making use of 
the inequalities ( |92] i and ( |94l i. 

B. Establishing ( |iS3l ) 

Fix n = 2, 3, . . . and consider = {K,p) with p in (0, 1), 
and positive integer such that iiT < n. 

Pick r = 1, 2, . . . , ri — 1. By Lemma [1 2. 11 we conclude that 



(95) 



Bn,r{0)Tn.l,---An,r < (1 - P) 



(96) 
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since < Un{9) < 1, and preconditioning arguments similar 
to the ones leading to dSTI l yield 



E* 



The event C„,,(^) depends only on r„i,...,r„r whereas 
E*^ is determined solely by r„^+i, . . . , r„ „. Thus, the 

event Cn,r{d) is independent of the rv (1 — p)^"-'' under the 
enforced assumptions, whence 



(97) 



Pick t arbitrary in (0, 1) and recall Lemma ri2. 21 A simple 
decomposition argument shows that 



E 



< E 



E* > (1 - t)rK 



n — 1 



E* < (1 - t)rK 



n — 1 

< e"(i-*)'-P^-^ + e 



Therefore, whenever r = 2, 3, . . . , [^J, we have 



E 



El 



since on that range we have 

n — r nil 1 
> — '- — > -. 

?i-l~n-l~2 

Now consider a scaling ; No — s> No x (0, 1) such that 
lini„_^ooPn = V* exists and (fTOl i holds for some c > t(j)*\ 
Replace d by 0„ in both ( |97] i and (|98] l according to this scaling 
and use the bound of Lemma [l2.4| in the resulting inequalities. 
Pick an integer i? > 2 (to be further specified shortly) and for 
n > 2{R + 1) note that 

n PK,.(0„)] 

r=R+l ^ ^ 

r=R+l ^ ^ 

LtJ 



< 



E 



r=R+l 

LtJ 



n- 1 



log n 



e 2 



•rpnKn 



E < 

r=_R+l 



^ 1 

n — 1 



loe; n 



by the same arguments as the ones leading to ( |92] i. Upon 
invoking the lower bound in ( |38] | we now conclude for all 
sufficiently large n > 2{R + 1) that 

LfJ 



E 

r=R+l 



< 



< 



E "^'^ 

r=_R+l 

LfJ 

+ E 

oo 

E 

r=_R+l 

oo 

+ E 



^ -1 

n — 1 



log 71 ) V 1°S" 



en ec 



n — 1 



logn I e 8 



— ^ -rCn log n 



n — 1 



log ri • n 



n — 1 



log n ■ n 



r=R.+ \ 

Furthermore, for all sufficiently large n > 2(i? + 1) it also the 
case that 

n 



n — 1 



1 f -i^c -tic \ 

logn -max In ^",n s ^" \ 



< 1 



(99) 



and the two infinite series converge. Let n*{R) denote any 
integer larger than 2(i? + 1) such that ( |99] l holds for all n > 
n*{R). On that range, by our earlier discussion we get 



LtJ 

E 



=i?+i 



'[AnA^n]] <e[ea 



n — 1 



R+l 



logn (...) 



with 



,(i?+i) 



1 - eCn:^\ogn -n 4 
„i-4c„(i?.+i) 

+ n^- 

1 - ec„;^logn ■ n s <=" 

Finally, let n go to infinity in this last expression: The 
desired conclusion ( [83] l follows whenever the conditions (1 — 
t)c{R+l) > 4 and c{R+l)i^ > 8 are satisfied. This can be 
achieved by taking R so that 

/ 4 

R + I > max 



This is always feasible for any given t in (0, 1) by taking R 
sufficiently large. ■ 



Appendix A 
A PROOF OF Proposition I1Q.1I 

The basis for deriving (ISST i lies in the observation that nodes 
1 and 2 are both isolated in H n G(n; 6) if and only if each 
edge in IEI(n; K) incident to one of these nodes is not present 
in G{n;p). Thus, Xn,i{(^) = XnA^) = 1 if ™d only if both 
sets of conditions 

Bi,{p) = if i]„,i n ^ 0, J e M-i 
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and 



hold. 



B2fe(p) = o if s„,2 n s„^fc ^ 0, fc e AA_2 



To formalize this observation, we introduce the random sets 
Nn^i{0) and Nn,2{0) defined by 

NnA^) {j = 3, . . . , n : j e r„,i V 1 £ r„,,} (100) 

and 

iVn,2(6') {fc = 3,...,n : fc G r„,2 V 2 G r„,fc}. (101) 

Thus, node j in Nn^i{d) is neither node 1 nor node 2, and is K- 
adjacent to node 1. Similarly, node k in Nn.2 (d) is neither node 
1 nor node 2, and is K-adjacent to node 2. Let Zn{0) denote 
the total number of edges in H(n; K) which are incident to 
either node 1 or node 2. It is plain that 

+ 1 [2 G r„4 V 1 G r„,2] (102) 

with the last term accounting for the possibility that nodes 
1 and 2 are K-adjacent. By conditioning on the rvs 
r„^i, . . . , r„^„, we readily conclude that 



E[x„4WXn,2W]-E[(l-p)^"('') 



(103) 



under the enforced independence of the collections of rvs 
{r„,i,...,r„,„} and <i<j<n}. 

To proceed we need to assess the various contributions to 
Z„{e): Using O we find 

n 

\Nnsm = ^i[jGr„.i viGr„,,] 

n n 

= ^i[jGr„a] + Ei[ier„,,] 



3=3 



3=3 



Gr„,i, 1 Gr„,, 



3=3 

n n 

= 51 1 [j e r„4] + 1 [.? ^ r„.i, 1 G r„,,] 

3=3 3=3 

= A'-l[2Gr„,i] 

n 

+ Y,i[j^Tns, 1 er„,,] (104) 

3=3 

where the last step used the fact |r„,i| = K. Similar 
arguments show that 



\Nn,2m = ^i[fcer„,2 v2Gr„,fc] 

fc=3 

= A--i[iGr„,2] 

n 

+ ^i[fc^r„,2, 2Gr„.fe] 

fc=3 



As a result, from the definition of Zn{9) we get 

z„{e) - 2 A- - 1 [2 G r„,i, 1 G r„,2] + z*{e) 



(105) 



(106) 



upon using ^ one more time, where 

n 

z^ie) ^i[j-0r„,i, iGr„ 



3=3 



+ ^i[.?^r„,2, 2Gr„j]. (107) 

3=3 

In order to evaluate the expression (1103) , we first compute 
the conditional expectation 



E 



(1-p) 



(108) 



From (IIO6I 1 we see that this quantity can be evaluated as the 
product of the two terms 



(1-p) 



2A'-(i[2er„,i4er„,2] 



and 



(109) 



(110) 



To evaluate this last conditional expectation, for each j = 
3, . . . , 71, we set 



Vn^j{e;S,T) 



with S and T subsets of M, each being of size K. It is 
straightforward to check that 

= 1 [j ^ S] 1 [j ^ T] E [(1 - p)i[i6r„,,]+i[2Gr„,] 



-l[j (^T]l[j eS]E 

-1 [j G S] 1 [.7 G T] . 



(l_p)i[ier„,,] 
(l-p)i[2er„,,] 



Then, with the notation introduced earlier in Section |IX] we 
can write 

Vn^,{9;S,T) 

= l[j ^S]l[j ^T]bn{9) 

+ (1 [j ^ S] 1 [j G T] + 1 [j ^ T] 1 [.? G S]) u^e) 

+1 [.? e S] 1 [j G T] . 

Next, the two rvs r„ 1 and r„ 2 being jointly independent 
of the rvs r„.3, . . . , r„.„, we find 



[ii-pr- 



rn,i7 rn,2 



(111) 



where the rvs Bn{9) and Un{d) are given by ( |56] l and ( |57| i, 
respectively. Therefore, since 



rri,i, rn.2 
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by a standard preconditioning argument, we get the expression 
dSST l upon writing ( 1108b as the product of the quantities (|1091 l 
and ( II 10b . and using ( llUb . ■ 



Appendix B 
A PROOF OF Lemma [T2TT] 

The defining conditions for Bn,r{&) lead to the representa- 
tion 

where we have set 

En.A^) {[k i r„,,] n ^ r„,fe]) u [B,k{p) = o] 

with i = 1, . . . , r and fc = r + 1, . . . , n. In terms of indicator 
functions, with the help of ([TJ this definition reads 

1 [EnMS)] 

= 1 [fc ^ r„.,] 1 [i i r„,fc] + (1 - B,u{v)) 
-i[fc^r„.,]i[i ^r„.fe](i-B,fe(p)) 
= (1 - + 1 [fc ^ r„,,] r„,fc] B,k{v)- 

Therefore, under the enforced independence assumptions, 

IP -B„^r(^') r„4, . . . , r„_„ 



Next, we find 

Gn,r{S\^ . . . , S'i 

= E 



n r 



= E 



■Q _p-)-Er=ii[fc^s',]i[»^r„.fc 

=r+l 



fc=r+l 

as we again use the enforced independence assumptions. Fix 
A; = r + 1 , ■ • • , and note that 



E 



n((i-p) 



.1=1 



E 



_z=l /c=rH-l 



where 



W{x\-p) ~\ - -p ^ ■px, a; G R. 
Since VF(.T,p) = (1 — p)^^'^ for .t = 0, 1, we obtain 

p[B„.,(0)|r„,i,...,r 

= E 



r n 



n n (1 ^) 

,i=l fc=r+l 



i-i[fc^r„,,]i[j^r„,fc] 



and it is now plain that 



r^i, 1 5 ■ • ■ ! -Tn.r 



= (i-py("-'-)G„,,(r„4,...,r„.,;0; 

where we have set 

Gn,r{Sl, ■ . ■ , Sr',d) 
r n 

.1=1 fc=r+l 

with Si, . . . ,Sr subsets of Af, each of size K. 



< 



r 

1[e[{i-p)-' 



l[k(^Si] 



(112) 



i=i 

where (II 12b follows from the negative association of the rvs 
- Use dHJ and note that 

^l^py^MS,>l^ z = l,...,r. 
Next we observe that for each i = 1, . . . , r, we have 



E 



(i-pY 

= (i-p)"'p[i ^r„.fe]+p[ier„,fc] 



A' 



a: 



n — 1 / n — 1 



1-p 



whence 



r 



Combining these observations readily yields 



k=r+l ^ ^ ^ 



Unjd) 
1-p 



We finally obtain 



^n, 1 ; ■ • ■ : rn,r 



Bn.rid) 

< (l_p)'-("-) 

and the desired conclusion ( |74] | follows. 



19 



Acknowledgment 
This work was supported by NSF Grant CCF-07290. 



References 



[1] 



I. F. Akyildiz, Y. Sankarsubramaniam, W. Su and E. Cayirci, "Wireless 
sensor networks: A survey," Computer Networks 38, pp. 393-422. 
N. P. Anthapadmanabhan and A. M. Makowski, "On the absence of 
isolated nodes in wireless ad-hoc networks with unreliable links - A 
curious gap," Proceedings of IEEE Infocom 2010, San Diego (CA), 
March 2010. 

S.R. Blackburn and S. Gerke, "Connectivity of the uniform random 
intersection graph," Discrete Mathematics 309 (2009), pp. 5130-5140. 
M. Bloznelis, J. Jaworski and K. Rybarczyk, "Component evolution in 
a secure wireless sensor network," Networks 53 (2009), pp. 19-26. 
B. Bollobas, Random Graphs, Second Edition, Cambridge Studies in 
Advanced Mathematics, Cambridge University Press, Cambridge (UK), 
2001. 

S. A. Qamtepe and B. Yener, "Key Distribution Mechanisms for Wireless 
Sensor Networks: a Survey," Technical Report TR-05-07, Computer 
Science Department, Rensselaer Polytechnic Institute, Troy (NY), Maixh 
2005. 

H. Chan, A. Perrig and D. Song, "Random key predistribution schemes 
for sensor networks," Proceedings of the 2003 IEEE Symposium on 
Research in Security and Privacy (SP 2003), Oakland (CA), May 2003, 
pp. 197-213. 

R. Di Pietro, L.V. Mancini, A. Mei, A. Panconesi and J. Radhakrish- 
nan, "Redoubtable sensor networks," ACM Transactions on Information 
Systems Security TISSEC 11 (2008), pp. 1-22. 

W. Du, J. Deng, Y.S. Han and P.K. Varshney, "A pairwise key pre- 
distribution scheme for wireless sensor networks," Proceedings of the 
10th ACM Conference on Computer and Communications Security 
(CCS 2003), Washington (DC), October 2003, pp. 42-51. 
D. Dubhashi and A. Panconesi, Concentration of Measure for the Anal- 
ysis of Randomized Algorithms, Cambridge University Press, New York 
(NY), 2009. 

L. Eschenauer and V.D. Gligor, "A key-management scheme for dis- 
tributed sensor networks," Proceedings of the ACM Conference on 
Computer and Communications Secuiity (CSS 2002), Washington (DC), 
November 2002, pp. 41-47. 

P. Gupta and P. R. Kumar, "Critical power for asymptotic connectivity 
in wireless networks. Chapter in Analysis, Control, Optimization and 
Applications: A Volume in Honor of W.H. Fleming, Edited by W.M. 
McEneany, G. Yin and Q. Zhang, Birkhauser, Boston (MA), 1998. 
J. Hwang and Y. Kim, "Revisiting random key pre-distribution schemes 
for wireless sensor networks," Proceedings of the Second ACM Work- 
shop on Security of Ad Hoc And Sensor Networks (SASN 2004), 
Washington (DC), October 2004. 

S. Janson, T. Luczak and A. Ruciriski, Random Graphs, Wiley- 
Interscience Series in Discrete Mathematics and Optimization, John 
Wiley & Sons, 2000. 
[15] K. Joag-Dev and F. Proschan, "Negative association of random variables, 
with applications," The Annals of Statistics 11 (1983), pp. 266-295 
G.E. Martin, Counting: The Art of Enumerative Combinatorics, Springer 
Veriag New York, 2001. 

A. Mei, A. Panconesi and J. Radhakrishnan, "Unassailable sensor net- 
works," Proceedings of the 4th International Conference on Security and 
Privacy in Communication Networks (SecureComm), Istanbul (Turkey), 
September 2008. 

K. Rybarczyk, "Diameter, connectivity and phase transition of the 
uniform random intersection graph," Submitted to Discrete Mathematics, 
July 2009. 

M.D. Penrose, Random Geometric Graphs, Oxford Studies in Probability 
S, Oxford University Press, New York (NY), 2003. 
A. Perrig, J. Stankovic and D. Wagner, "Security in wireless sensor 
networks," Communications of the ACM 47 (2004), pp. 53-57. 
D.-M. Sun and B. He, "Review of key management mechanisms in 
wireless sensor networks," Acta AutomaU'ca Sinica 12 (2006), pp. 900- 
906. 

O. Yagan and A.M. Makowski, "On the random graph induced by a 
random key predistribution scheme under full visibility," Proceedings of 
the IEEE International Symposium on Information Theory (ISIT 2008), 
Toronto (ON), June 2008. 



[2] 



[3] 
[4] 



[5] 



[6] 



[7] 



[8] 



[9] 



[10] 



[11] 



[12] 



[13] 



[14] 



[16] 
[17] 



[18] 



[19] 



[20] 



[21] 



[22] 



[23] O. Yagan and A. M. Makowski, "Connectivity results for random 
key graphs," Proceedings of the IEEE International Symposium on 
Information Theory (ISIT 2009), Seoul (Korea), June 2009. 

[24] O. Yagan and A.M. Makowski, "Zero-one laws for connectivity 
in random key graphs," Available online at arXiv:0908.3644 
[math. CO], August 2009. Earlier draft available online at 
http://hdl.handle.net/1903/8716 January 2009. 

[25] O. Yagan and A. M. Makowski, "On random graphs associated with 
a pairwise key distribution scheme for wireless sensor networks (Ex- 
tended version)," submitted for inclusion in the program of IEEE 
Infocom 2011, Shanghai (PRC), April 2011. Available online at 
http://hdl.handle.net/l 903/ 1 060 1 

[26] O. Yagan and A. M. Makowski, "On the gradual deployment of random 
pairwise key distribution schemes," submitted for inclusion in the 
program of IEEE Infocom 2011, Shanghai (PRC), April 2011. Available 
online at http://hdl.handle.net/1903/10604 

[27] O. Yagan and A. M. Makowski, "Designing securely connected wireless 
sensor networks in the presence of unreliable hnks," submitted for 
inclusion in the program of ICC 2011, Tokyo (Japan), June 2011. 

[28] C.W. Yi, PJ. Wan, K.W. Lin and C.H. Huang, "Asymptotic distribution 
of the number of isolated nodes in wireless ad hoc networks with 
unreliable nodes and links," Proceedings of IEEE Globecom 2006, San 
Francisco (CA), November 2006. 



